Numerical method for finding decoherence-free subspaces and its applications 
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In this work, inspired by the study of semidefinite programming for block-diagonahzing matrix *- 
algebras, we propose an algorithm that can find the algebraic structure of decoherence-free subspaces 
(DFS's) for a given noisy quantum channel. We prove that this algorithm will work for all cases 
with probability one, and it is more efficient than the algorithm proposed by Holbrook, Kribs, 
and Lafiamme [Quant. Inf. Proc. 80, 381 (2003)]. In fact, our results reveal that this previous 
algorithm only works for special cases. As an application, we discuss how this method can be applied 
to increase the efficiency of an optimization procedure for finding an approximate DPS. 
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I. INTRODUCTION 

Decohercnce and other noises cause errors in quan- 
tum information processing. Several methods have been 
proposed to significantly reduce these errors [l|, such 
as quantum error correcting codes 2|-|8|, decoherence- 
free/noiseless subsysten is ( DFS's) |9l-[l9|. and dynami- 
cal decoupling (DD) j2(]| - |2a |. Hovifever, dynamical decou- 
pling controls often do not satisfy the strict experimental 
requirements and quantum error correcting codes, while 
possibly protecting against any error, require a large re- 
source overhead. Decoherence-free subspaces, or noise- 
less subsystems, can reduce overhead since they do not 
require detecting and correcting errors but are difficult 
to identify and use. 

In terms of the (Krauss) operator-sum representation, 
the noisy quantum evolution can be fully characterized 
by a set of operators that generate a noise algebra A. 
The structure of a DFS can be recovered by studying 
the algebraic structure of A, or its commutant algebra 
A' , both of which are special examples of a so-called C*- 
algebra, or matrix *-algebra [l^l- Although DFS's can 
be obtained analytically for certain noisy systems, this 
is not possible in general. In [28| a numerical algorithm 
is proposed in which the noise algebra A is used to cal- 
culate the commutant algebra A! and decompose it into 
the algebraic form which gives the structure of all DFS's 
if any exist. However, we find that this algorithm is in- 
complete although useful for the special solutions of the 
basis of A' as chosen in [l^ . In practice, in most of the 
cases, computers will pick up other solutions in which 
their algorithm does not give the complete decomposi- 
tion. Because of this we were motivated to find a new 
algorithm that provides a general algorithm which can 
be explicitly implemented to give the required decompo- 
sition. 

In fact, since A and A' are special examples of a matrix 
^-algebra, we can try to solve the more general problem 
of how to decompose an arbitrary matrix *-algebra. We 
find that this problem is equivalent to that of simultane- 



ously block-diagonalizing a matrix *-algebra, and this has 
been well-studied in research on semidefinite program- 
ming [29l. Isoj. In particular, in j29| . a numerical method 
was proposed to find the finest block-diagonalization of 
the algebra generated by real symmetric matrices. Their 
method consists of two steps. First, decompose the al- 
gebra into simple components, and second, decompose 
each simple component into irreducible components (the 
details of this procedure will be clarified in the follow- 
ing). In this work, we will show that such a two-step 
algorithm can also be applied to the algebra generated 
by Hermitian matrices, such as A and A' . We will also 
give analysis and proofs to show the validity of this algo- 
rithm. As applications, we apply our new algorithm to 
the collective- noise model in [28|, and compare the nu- 
merical results with the algorithm proposed there. We 
find that our improved algorithm is not only valid in 
general, but is also more efficient: it requires fewer con- 
ditional loops, and requires only A or AJ alone, rather 
than both [28[. 

The paper is organized as follows: in Section |TT1 we in- 
troduce the matrix *-algebra A generated by noise oper- 
ators for a given quantum channel, and the Wedderburn 
decomposition for a DFS. In Section HIIl we present an al- 
gorithm to transform A into the Wedderburn form using 
two steps. Finally in Section HVl for the collective-noise 
model, we numerically implement and compare our algo- 
rithm with the one proposed in [2^ . We will also briefly 
discuss how our algorithm can be used to find a good 
initial point for the optimization process in searching for 
an approximate DFS. 

II. ALGEBRAIC STRUCTURE OF A NOISE 
ALGEBRA 

A. Noise Algebra for a Noisy Quantum Channel 

Let p be the density operator of an n-dimensional 
quantum system with Hilbert space %. In real physical 
systems, the evolution of p suffers from noise due to its in- 
teraction with the environment. Such noisy evolution can 
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be represented as a quantum channel £: p ^ ^{p)i where 
£ is characterized by a set of operators {Aj.}, j = 1, . . . ,p: 



fc=l 



fc=2 



This is often referred to as the Kraus operator-sum rep- 
resentation (3ll |. where {A^} can include a Hamiltonian 
as well as irreversible coupling to a Markovian bath. In 
the following, we will assume the Hamiltonian H = 
and only focus on the noise effect on p. In this case the 
{Ak} contain information purely about the noise, and are 
known as noise operators. For many channels, we have 
Ai = y/Po^, where po represents the probability that no 
error occurs. The noisy channel is referred to as uni- 
tal if £{!) = Efc^fc4 = I- Define the noise algebra A 
to be the C*-algebra, or the matrix *-algebra generated 
by {Afe}. The definition of a matrix *-algebra simply 
implies that A is closed under matrix summation, multi- 
plication and f-operation. The reason why we introduce 
the concept of a matrix *-algebra is that it can be de- 
composed into a nice algebraic structure, with details in 
the following. 



B. Wedderburn Decomposition for a DFS 



Here Afi = Aim Im; is an algebra different from 
A^n, © • ■ ■ ffi A^n, • Applying Theorem [T] to the conju- 
gates A and A' , we can find some unitary U such that: 



UUU = = 07W„, ® (2a) 

i i 

e i 

C/U'C/ = 0A/;' = 0I„,0M™.. (2b) 



Mathematically, each Afi, i = 1, • • • corresponds to a 
simple component of A, while the subblock Mi at each 
diagonal position corresponds to an irreducible compo- 
nent. 

Assume that there exists some > 1, and call this 
nik- We can encode an arbitrary mfc-dimensional state p 
into p = I„j^ (g) p © Ores G A' such that £{p) = p, where 
Ores represents the zero density operator on the rest of 
the Hilbert space with respect to Op. Hence, if wc find 
the Wedderburn decomposition for Aor A' , then each Afi 
with mi > 1 corresponds to a decoherence-free subsystem 
(which reduces to a decoherence-free subspace if = 1). 
Moreover, since A and A' obey the commutant relation 
given in Eq. ([2]) , we do not need both the decompositions 
for A and A'; one will suffice. 



For a general p, £{p) 7^ p, and thus the quantum infor- 
mation stored in p will not be preserved by the noise. It 
may be possible, however, find a subspace or subsystem 
in some space Hi E H such that for p £ Hi, £{p) = p. If 
so "Hi is called a decoherence free subspace or subsystem 
(DFS). Notice that for a unital channel £, if [p, A^] = 0, 
for all k, then £{p) = p. Hence to locate a DFS it is 
enough to study the commutant of A, which is defined 
to be 

A' ^{B\[B,A] =:0, AeA}, 

and is also a matrix ^-al gebr a. Applying the 
Wedderburn-Artin theorem [27ll3^ to a special case, it 
can be shown that every matrix *-algebra with an iden- 
tity has the following fundamental structure decomposi- 
tion Hi,!!!]: 

Theorem 1. (Wedderburn decomposition) Let A C 
C"^" be a matrix ^-algebra with an identity. Then there 
exists a unitary transformation U such that U^AU has a 
block- diagonal structure: 

U^AU = diag(A/'i, A/'a, ■ ■ ■ Mt) 

where each Afi corresponds to a simple subalgebra com- 
ponent. Moreover, Afi has the following block- diagonal 
structure: 

M = {diag(M,, • • • , Mi), e AJ„J = AA^ ® Im. 

(1) 

where AAm denotes the Ui x Ui matrix ^-algebra over the 
complex field C. 



III. NUMERICAL ALGORITHM TO OBTAIN 
THE WEDDERBURN DECOMPOSITION 

To find the Wedderburn decomposition for a quantum 
channel given by a group of noise operators {^/c}, it is 
sufficient to find the unitary transform U such that A 
and A' are simultaneously block-diagonalized into the 
decomposition in Eq. An algorithm to do this for 
real symmetric Ak is given in [29j . Here we construct an 
equivalent algorithm that we prove works for Hermitian 
Ak- This is sufficient for our purposes, because while the 
noise operators Ak need not be Hermitian, we can always 
replace a non-Hcrmitian operator Aj with the two Her- 
mitian operators A^^' = Aj + A] and A^^ = iiA,-. - A\) 

J J J A 

and still have a generating set for the algebra A. For 
simplicity we simply assume that all the Ak are Hermi- 
tian and form a basis for A. The advantage of choosing 
an Hermitian basis will be clear in the following analysis. 
The algorithm breaks into two steps: 

Algorithm 1. (Wedderburn decomposition) Let A C 
C"^" be a matrix ^-algebra with an identity, and A a 
"generic" element of A ("generic" is defined below). 
Step 1: Find the unitary transform V such that 

yUy = diag(Ci,C2,--- ,Cf), (3) 

where each Ci corresponds to some representation of the 
simple component Afi in (0^. 

Step 2: Find the local unitary transform P such that 
within each Ci, P^V'' AVP is equal to Afi — Aim ® Imi • 
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Then U = VP is the required unitary transform for the 
Wedderhurn decomposition. 

To implement the two steps above one picks a single 
operator A ^ A, and diagonalizes A to find the required 
decompositions. Due to the decomposition in Eq. ([2]), we 
know that there exists a unitary transformation V such 
that for Ae A, 

I 

y\AV = 0(I„. ® A), (4) 

i 

where the Di are diagonal matrices whose elements are 
the eigenvalues of A. To obtain the spaces spanned by 
the simple algebras from this eigenvalue decomposition, 
we need to pick an A such that the Di,. . . ,Di do not 
share any eigenvalues, and the eigenvalues in each Di are 
distinct. Note that A will have this property if it has 
the maximum possible number of distinct eigenvalues. It 
can be shown [2^ that the set of operators that have 
this maximum number of distinct eigenvalues is topolog- 
ically dense in A, and so we will refer to such A as being 
generic. If we randomly choose A from A using a suitable 
measure, it will be generic with probability 1. A simple 
way to generate a generic A is to choose a random vector 
a. = {ai, . . . , a-k) and generate A ~ '^j^j- 

After picking a generic A, we diagonalize A and obtain 
the distinct eigenvalues, Xj, and their multiplicities, kj, 
j = 1, . . . ,q. We then group these eigenvalues and write 
down the eigenspace decomposition according to their 
multiplicities in a non-decreasing order: 

V^AV = diag(AiIfc, , Aalfe, , • • • , Xqlk, ) (5) 

with V = (Vi, V2, . . . , Vq) where each Vj is composed 
of the eigenvectors corresponding to the eigenspace of 
Xj. We can further define a new division of V: V = 
J<i-')) where each ifW is the union of 
all eigenspaces Vj with the same multiplicity p,.: 

K^^^ = V, 

kj=Pr 

where pi < P2 < ■ ■ ■ < Ps are the distinct multiplicities 
of the eigenvalue a/s. 

By Theorem [TJ each Vj will lie in some simple compo- 
nent A/i, I = 1, . . . , ^. Due to the form of Ni = Mm ®^mi , 
we immediately know that only V^'s within the same 
can belong to the same A/i, and Vj's in differ- 
ent K^'^'^^s must belong to different A/i's. Hence, each 
_ftr'^^ must either be some A/i, or a direct sum of a few 
A/i's, in which case can be further decomposed. 

In either case A is block-diagonalized over the division 

There is a simple method to check whether if'^'"^ can 
be further decomposed: wc choose another randomly 
generated A — Pj^j ^ Since a and /3 are 

independent, with probability 1, A and A will gener- 
ate the whole algebra A. As we have pointed out. 



both A and A are block-diagonalized over the division 
,K^'')). If there exist Vj and Vy within 

some K^''''> such that VjAVf = 0, then A can be further 
block-diagonalized on A''*"^ over the division between Vj 
and Vj'. In this way, by checking the value of Vj AVj' 
between all different j and j' on each K^'^\ we can iden- 
tify the structure of each Ai in each K'^^\ and finally 
make both A and A simultaneously block-diagonalized 
over (BiMi. Since A and A will generate A with proba- 
bility 1, we can claim that the whole algebra A has been 
simultaneously block-diagonalized over OiAfi. However, 
we should notice that within each sub-block Ai , A may 
not be the same as A^„. ® I^^, but some representation 
of it, so we will instead denote the sub-block by C,. Thus 
we have obtained a V that transforms A into the form 
of Eq. ([3]). In particular this V has already transformed 
the generic A into the Wedderburn form: 

e 

V^AV = ^{D,<i^I^^), (6) 

i 

where Di is a diagonal matrix with all distinct eigenvalues 
of A on each Ci . 

Now we note that V''AV is usually not in the form of 
A/i on Ci. In the next step, we are looking for a further 
unitary transform that leaves V'' AV invariant but trans- 
forms V^AV into the form of A/i on each Ci. Without 
loss of generality, let us focus on a simple component Ci 
which is composed of a few eigenspaces Vj of A: 

C, = Vl'^ ® V^^'^ (B---®Vj;i} 

If according to the division ©iCi, we define a local unitary 
transform P to be: 

P = diag(p(i\p(i),--- ,P(^)) (7) 

P«^diag(P«,P«,... ,P«) (8) 

where Pj'*'' is a unitary matrix on the subspace vj'^'^ , then 

such P will leave V''AV invariant. Moreover, the follow- 
ing result is proved in Proposition 3.7 in [29j : 

Theorem 2. For A and V satisfying (0), there ex- 
ists a local unitary transform P as in ^ such that 
PW^AVP ^®,{Mn.®l^^)- 

Hence, it is possible to construct a local unitary Q 
(which may not be equal to P) in the form of Eq. ([7]) 
such that Q^V'^ AVQ is in the Wedderburn form. Before 
we design the required Q, we would like to find out what 
the matrix of V^AV looks like on C; after the transform 
V. Notice that since Theorem [5] only claims the existence 
of such P, the local transform Q we finally construct may 
look either the same as, or different from P. 

A matrix is called a scalar matrix if it is equal to a 
scalar times an identity matrix. On each Ci, as a corollary 
of Theorem ^ we have 

A,,, ^ (V^AV)^, = vj'»AvjP = h,P^'^P^P' (9) 
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For j = j\ Aj,ji is equal to kjjlrm (that is, the di- 
agonal sub-blocks of (V^^Ay)'*) are already in scalar- 
matrix form). For j ^ j' the off-diagonal sub-blocks 
Ajji may or may not be in this form. Our next goal is 
to find a local transform Q in the form of (0 such that 
((5(*)t(yt^y)(i)g(i))^.^., are in scalar-matrix form for all 
j and f. 

(i) 

On each d, we can sequentially construct each Qj 

in Q(*) EE diag(Q^'\ Q^*\ • • • , Qt\). First, choose Q^^ = 
Im. . Then for j > 2 define 



where qj is the first row of Q^*^ . We now prove that this 
Q'*^ is the required unitary transform. 

Theorem 3. Q^*^ as defined in U0\} are unitary matri- 
ces, and Q^*^ = diag((5^*'', \ ' ' ' unitary 
transform such that Q^*'' {V^ AV)^^''j,Q^^) are scalar ma- 
trices. 

(i) 

Proof. To show Qj is a unitary matrix, it is sufficient to 

show Qj Qj is in scalar-matrix form. For j > 2, from 
([9|), we have: 

Hence, after normalization, Q^'^ becomes a unitary ma- 
trix. In addition, 

Qf{V^AV)'^,Q'P = Qf^'fc,-,'P«P«^Q« 



(10a) 
(10b) 



,(Otp('i 

i' j' 



and so all sub-blocks are in the scalar-matrix form. 



□ 



Numerically, following Eq. (|TOl). we can construct the 
local unitary transform Q in the form of Eq. ([7]) that leave 
A invariant but transforms A into the form ©i(M„; (g) 
ImJ- Since A is generated by A and A, we can claim 
that the whole algebra A is in the Wedderburn form after 
the unitary transform U = VQ. We summarize the full 
algorithm in Table H] 

When implementing the algorithm for a given noisy 
channel, we can calculate the Wedderburn form for ei- 
ther A or A' , depending on which one is easier to derive. 
Notice that for special cases when Mi = A^a_, (g) Ii, or 
Ni = Ik, after Step 1 in Algorithm 1, Ci will already be 
the same as Afi. For such cases, there is no need to im- 
plement Step 2, and we can simply choose Q*-*' = I on 



step 1: (a) from A, pick two generic matrices A and A 

(b) diagonalize A and A to get V as in Eq. ((5]) 

(c) find the structure of A/i in K^''\ getting Eq. (|6} 
Step 2: (d) build tlie local transform Q using Eq. ((TOD 

(e) U = VQ is the required unitary in Eq. ^ 

TABLE I. Algorithm to find U in the decomposition Eq. ©. 



IV. APPLICATIONS 

A. Finding the DFS Structure of a Channel 

The primary application for the above algorithm is de- 
riving the DFS structure for a given noisy quantum chan- 
nel, and finding the corresponding unitary transform U 
in Eq. First of all, let's reinvestigate the collective 
noise model calculated in 



28j . For a system with Uq 



qubits, we say a quantum channel £ is under collective 
noise if 

£{p) = A^pAl + AypAl + A.pAl 
^fc = -^e*'^", k = x,y,z, 

where 



Sx — ^ ] Sy — ^ ^ Yi, Sz — ^ ] Zi 

i—l i—1 i=l 

are sums of local Pauli operators on each qubit. For such 
a noisy channel, we can define the algebra generated by 
the noise operators by 

A EE span{Ax, Ay, A^} = span{S'^, Sy, S^}, 

= span{So,, Sy} = spanj^j^, S^} = spanlS*:!;, S^}, 

We would like to find the Wedderburn decomposition 
of A or A' in the form of Eq. Notice that since 

the collective noise channel is a special type of noisy 
channel, we can actually derive the the fundamental de- 
composition by using Young diagrams for the addition 



of angular momentum for any value of Uq [34l |. For ex- 
ample, for Uq = 3, A ^ {M2 (8) II2) ® Mr, for = 4, 
A = I2© (A^s^Is)® A^s- However, theory does not give 
a specific basis for the operators A, and must identify one 
numerically. In the following, we shall apply both the al- 
gorithm suggested in and the above Algorithm 1 to 
the collective noise channel and compare the numerical 
results. 

In the algorithm suggested by [l^, we need to first 
calculate each Bj in A' = spanjPi, . . . ,Br}, where {Bj} 
form a basis for A' . Next, based on the operators {Bj}, 
j = ,r, we find a group of so-called minimal- 

reducing projectors Pj, and then block-diagonalizc A' 
into the form diag(Ci,C2, ■ ■ • ,Ci), where Ci is a repre- 
sentation of Afi = Aim <8) Inii- Then after reshuffiing 
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the order of basis vectors, all the diagonal sub-blocks 
are transformed into the scalar-matrix form. It is then 
claimed in that the whole algebra A' is in the form 
of Eq. ©. 

If we compare the algorithm in [2^ with the algorithm 
in Table U we find that the former achieves Step 1, but 
Step 2 is missing. Step 2 turns out to be necessary for 
most cases, since the solution {Bj} as the set of basis 
A' is not unique. It is true that for the {Bj} chosen 
in [H] , Step 1 and reshuffling of basis vectors are enough 
to transform A' into the form of ([2]), but such choice of 
{Bj} is very special. A different solution for the {Bj} 
will be obtained it is derived numerically by solving the 
system of linear equations [Bj, ^fc] = 0, A; = 1, • • • ,p. 

As an example of the necessity of step 2, consider 
for riq = 4, and A' = M2 ® {M3 «) I3) ® Is- Follow- 
ing the algorithm in [2^ and using Matlab, we derive a 
group of 14 orthonormal basis matrices in A' {{Bj}, j = 
1, • • • , 14) which are different from those in Then we 
find the corresponding minimal- reducing projectors P,„, 



m = 1, • • ■ ,6, in which rank(Pi) = rank(P2) = 1 corre- 
sponding to the Afi ~ M.2 subspace, rank(P6) = 5 cor- 
responding to the As — I5 subspace, and rank(Pi) = 3, 
i = 3, 4, 5, corresponding to the J\f2 = M3 ® I3 subspace. 
Hence, {Pm} induces a unitary transform V such that 
V^A'V is in the form diag(Ci, C2, C3), where Ci = TVi, 
C3 = As, and C2 is some representation of N2- After 
reshuffling the basis vectors we can make the three di- 
agonal sub-blocks of V'^ BjV on C2 in the scalar- matrix 
form. 

Specifically, if we still denote BiV as Pi, then after 
reshuffling, on C2 we have 

/-O.I67I3 Pl,2 Pl,3 \ 

Pi = PI,2 O.233I3 P2,3 

V 51,3 ^2,3 0.308l3,y 

where 



-0.180 -h0.089z -0.120- 0.242i -0.042 0.063i> 
Pi2=| 0.097 + 0.169« 0.103 -O.OOH -0.259 0.059i 
0.193 -O.OeOi -0.100- 0.159i -0.051 - 0.200i; 

'-0.074 + 0.136i 0.055 + 0.047i 0.072 + 0.055i 
Pi3= I -0.003 + 0.087i -0.094 + 0.080i -0.074 - 0.095i 
-0.066 + 0.039i 0.087- 0.09601 -0.114 - 0.041iy 

'-0.030- 0.078i -0.020- 0.068i -0.067 + 0.180^^ 
^2,3 = I -0.118-h0.099j 0.139 -hO.OSOi -0.037 0.045i 
^-0.133- 0.023i -0.016- 0.147i -0.009 - 0.093iy 

Therefore, we see that for a general solution {Bj} for A', such as the solution derived from the matlab routine, the 
algorithm suggested by [l^ fails to give the Weddcrburn decomposition, although it does give the the correct form 
for the particularly chosen {Bj} in [28]. Hence, in practice, the algorithm in [28[ is sometimes insufficient, which 
motivated us to develop the modified algorithm presented here. Next, wc would like to continue with this example 
for Uq = 4, and following our algorithm to find the local unitary transform Q = diag((5i, (52, <53) such that all the 
sub-blocks of Q'^BjQ are in the scalar matrix form. Strictly speaking, wc should use the random combination method 
in the last section to pick up a generic P for Step 2. However, for this particular example, the above Pi is already a 
generic matrix, so we will instead based on Pi to construct Q. 

Define Qi ~ I3, and according to Step 2 in TableU we define Qk = B^l and Qk = l/||<Zfc||Qfc, fc = 2, 3, where qk is 
the first row of Qfc. Then we have 

-O.I67I3 O.345I3 O.I93II3 

g^PiQ= ( O.345I3 O.233I3 -0.033 - 0.219il3 

O.193I3 -0.033 0.219«l3 O.3O8I3 



We can double-check the form of Q^BjQ for other Bj 
and we will find that all Bj are transformed in the form 
of A^3 0/3. Hence, for this particular set of {Bj}, we 
have explicitly constructed the unitary matrix U = VQ 
that transforms A' into the Weddcrburn decomposition, 
which is not accessible from the algorithm in l2g|. Now 



we know how to encode an arbitrary three-level quantum 
state p into the DFS J\f2- In the current computational 
basis, the encoded density operator p should take the 
form: 

p = c/(oe (p® I3) ® o)pt 

We note finally that: i) in the above implementation of 
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our algorithm, we have skipped Step 1 since the algorithm 
in has already transformed all Bj into [3] ii) We have 
applied Step 2 to Bi instead of to a random combination 
of the {Bj}, so as to make it easy for comparison. In 
practice, we do not really have to use the two generic 
operators as we suggest. They are introduced primarily 
to guarantee the validity of the algorithm. 

B. Searching for an Approximate DFS 

Although for every noisy channel there exists a decom- 
position as in ([2]), not all channels have a DFS that is 
useful for protecting quantum information. In fact, the 
noisy channels with a useful DFS constitute only a very 
small set. For many channels the algebraic decomposi- 
tion ([2]) looks like the following: 

i 

where ki ^ kj, for i j, which means that all the A^„.'s 
are 1-dimensional and so cannot store quantum informa- 
tion. When this happens, we would like to ask an alter- 
native question: does there exist a subsystem on which 
the noise, even if not zero, is significantly reduced? This 
is the concept of an approximate DFS (ADFS). 

It is not easy to characterize an ADFS by algebraic 
conditions, as we have done for a perfect DFS. Rather, 
an ADFS should be formulated as an optimal solution 
such that the noise on the system is reduced as much as 
possible. Hence, it is possible to obtain an ADFS numeri- 
cally by solving the corresponding optimization problem. 
Since there is more than one way to quantify the effects of 
noise, there is more than one way to define the function 
to be minimized in the optimization. For the purposes 
of our analysis in what follows, we simply assume that 
one such function has been chosen, and denote it by J. 
Furthermore, one must specify when the noise is "small 
enough" to be helpful as an ADFS. 

The problem of finding an ADFS involves searching for 
the optimal unitary matrix U that transforms the origi- 
nal basis into a new basis, such that a state pi, encoded 
in p = ® I„2 ©Ores, experiences the least noise under 

the noise operators {Ak}, where Ak = UAkW- That is, 
we want to minimize J{U), where U varies over the uni- 
tary group. Numerically, we can ap ply the BFGS quasi 
Newton method for the optimization |35| . Broadly speak- 
ing this involves i) choosing an initial point U — U^'^'^ in 
the unitary matrix space; ii) calculating the value and 
the gradient of the objective function J*^"^ = J[?7(°)]; iii) 
using the value and the gradient to implicitly derive the 
Hessian information and use them all to get a new 

C/(i) 

such that J[t/(^-'] < J[J7'-°^]. Repeating steps ii) and iii), 
we obtain a sequence of {U^'^^} with a limit corresponding 
to a local minimum of J. The limiting unitary transform 
U is what we are looking for. 

Due to the existence of many local minima, different 
initial choices of U^'^'^ may result in different values of the 



N 


1 


2 


3 


4 


5 


6 


J min 


0.2195 


0.2078 


0.2088 


0.0123 


0.4788 


0.0125 



TABLE II. The final value of Jmin obtained from numerical 
optimization, using 6 different random starting points enu- 
merated by A'^. 

optimized J. In particular, as the dimension gets larger 
we may need to run the optimization many times be- 
fore we obtain a value of J close to the global minimum. 
Hence a wise choice of U^^"^ can be very important in per- 
forming the numerical optimization. One important re- 
sult in optimization theory is that any gradient-based al- 
gorithm only guarantees that the iteration sequence will 
converge to some local minimum; however, if the initial 
point of optimization iteration is very close to the global 
minimum, then the iteration sequence will converge to 
the global minimum. On the other hand, if our noisy 
quantum channel £ can be considered as a perturbation 
of another channel £' that has a perfect DFS, then the 
ADFS of £ should be pretty close to the DFS oi£. Hence, 
we can apply the DFS algorithm in Scction lHll to first cal- 
culate the unitary matrix Uq for the DFS of £' , and then 
run the optimization for ADFS, with U^^^ = Uq- In that 
way, we will be able to derive the ADFS more efficiently. 

Specifically, let us take the collective noise model as an 
example, but this time add a perturbation to the original 
noise operators: 

Ax = V^Ax, Ay = Ay, Az = Az, 

where we have defined a perturbation unitary matrix 
that is sufficiently close to the identity: — I|| < e. 
Under the new noise operators Ax^y^z, we apply the al- 
gorithm in Section IIIII and find that there is no useful 
DFS. Next we try to run optimization to this model in 
searching for ADFS. First we do the optimization using 
random initial point. For example, choosing e = 1 and 
Uq = 4, we run the optimization routine 6 times, starting 
from different initial U^^\ and record the final minimized 
Jmin in Table [Hi We see that among the six different 
runs, only two of them have obtained Jmin < 0.0125. 
Hence we cannot guarantee that we have the best min- 
imized J from a single run of the optimization process 
from an arbitrary random initial {U^^^}. However, if we 
instead choose J7'-°^ = Uq, where Uq is the unitary matrix 
in the Wedderburn decomposition for the perfect DFS of 
{Ax.y.z}, then the optimization generates the minimized 
Jmin ~ 0.0123. For other values of Ug we find similar 
results. Thus our DFS-finding algorithm helps in finding 
good initial points for the optimization of ADFS search- 
ing. 

In addition, many local minima may also result in op- 
tions for our ADFS implementation. Whereas our algo- 
rithm gives the dimensionally optimal DFS, there may 
be several ADFSs and the "best" one may not be the 
dimensionally optimal one. The "best" might be the one 
which has robust, experimentally available controls, or 
one that has the lowest error rate per unit time. 
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